fj APPLICATION 

1! FOR 

g UNITED STATES LETTERS PATENT 

m Be it known that we, Alexander Berk, residing at 8 Canoe River Road, Sharon, MA 

5 i|| 02067, Prabhat K. Acharya, residing at 16 Upland Rd, Burlington, MA 01803, Lawrence S. 

Q Bernstein, residing at 6 Turning Mill Road, Lexington, MA 02420, Gail P. Anderson, residing at 

875 7 th Street, Boulder, CO 80302, Paul Lewis residing at 6200 Wilson Blvd, Apt. 1003, Falls 

Church, VA 22044, James H. Chetwynd, residing at 57 Norval Avenue, Stoneham, MA 02180 

and Michael L. Hoke, residing at 65 Long Hill Dr. Leominster MA 01453, and being a citizens 

10 of the United States of America, have invented a certain new and useful 

BAND MODEL METHOD FOR MODELING ATMOSPHERIC PROPAGATION AT 
ARBITRARILY FINE SPECTRAL RESOLUTION 

of which the following is a specification: 
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Applicant: Berk et al. 

For: Band Model Method for Modeling Atmospheric Propagation at Arbitrarily Fine 

Spectral Resolution 

GOVERNMENT RIGHTS 

5 This invention was made with Government support under Contract No. F19628-98-C- 

0050 awarded by the Department of the Air Force. The Government has certain rights in this 
invention. 

FIELD OF THE INVENTION 

This invention relates to a radiative transport band model for application to high spectral 
1 0 resolution simulations. 

BACKGROUND OF THE INVENTION 

Jj MODTRAN4 is the U.S. Air Force (USAF) standard moderate (2 cm -1 ) or 

Ill broader/coarser spectral resolution radiative transport (RT) model for wavelengths extending 

15 5 from the thermal InfraRed (IR) through the visible and into the ultraviolet (0.2 to 10,000.0 jjm). 
^ [See: A. Berk, G.P. Anderson, P.K. Acharya, L.S. Bernstein, J.H. Chetwynd, M.W. Matthew, 
% E.P. Shettle and S.M. Adler-Golden, "MODTRAN4 User's Manual Air Force Research 
jfj Laboratory Report, June 1999, and see: A. Berk, L.S. Bernstein, G.P. Anderson, P.K. Acharya, 
H DC. Robertson, J.H. Chetwynd and S.M. Adler-Golden, cc MODTRAN Cloud and Multiple 

20 Scattering Upgrades with Application to AVIRIS," Remote Sens. Environ. 65, pp. 367-375, 
1998, the disclosures of which are incorporated herein by reference.] 

The MODTRAN4 1 cm* 1 statistical band model (from which 2 cm" 1 spectral resolution 
results are obtained) provides a fast alternative (100-fold increase in speed) to the USAF more 
accurate line-by-line (LBL) radiative transport models, FASCODE and FASE. [See: S. A. 
25 Clough, F X. Kneizys, G.P. Anderson, E.P. Shettle, J.H. Chetwynd, L.W. Abreu, and L A. Hall, 
"FASCODE3 Spectral Simulation," Proceedings of the International Radiation Symposium^ 
Lenoble and Geleyn, Deepak Publishing, 1988, and see: H E. Snell, J.-L. Moncet, G.P. 
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Anderson, J.H. Chetwynd, S. Miller, and J. Wang, 'TASCODE for the Environment (FASE)", 
Proceedings of Atmospheric Propagation and Remote Sensing IV, SPIE, 2471, pp. 88-95, 1995.] 
FASCODE and FASE are both based on a 'first principles' physical equations, expanding the 
optical depth terms, based on spectroscopic constants and line shape, with very high accuracy. 
Comparisons between MODTRAN4 and FASE spectral transmittances and radiances show 
agreement to within a few percent or better in the thermal IR. MODTRAN4 includes flux and 
atmosphere-scattered solar calculations, essential components in analysis of near-IR, visible and 
ultraviolet spectral region data that are not readily generated by LBL models. 

MODTRAN4 and its predecessors have been used extensively over the last quarter 
century in the design and analysis of broadband, multiband, and short-wave IR / Visible 
hyperspectral imaging sensors. However, conventional interferometers and many state-of-the-art 
sensors working in the long- and mid-wave IR operate at higher spectral resolution than 
MODTRAN4 provides. 

Narrowing the band model spectral resolution changes the fundamental character of the 
band model. The half-width of molecular transitions near sea level average about 0.08 cm" 1 . As 
illustrated in Figure 1, the 1.0 cm" 1 band model calculates the absorption of atomic and molecular 
lines whose line center regions lie almost entirely within the spectral bin. At the finer spectral 
resolution, a much larger fraction of any atomic or molecular line falls outside of the spectral bin 
containing the line. Determination of the new band model has therefore required improved 
treatment of both line tail and line center absorption. Line tail absorption is modeled closer to 
line centers (as defined by a compilation of spectroscopic data), and the finite-bin single-line 
equivalent width used to calculate line center absorption is no longer simply a small perturbation 
of the total single line equivalent width. 
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The line center absorption within a spectral bin is generally defined as the in-band 
absorption from all atomic and molecular transitions centered in that bin, Figure 2. LBL models 
calculate this in-band absorption by explicitly determining the spectral absorption of each line on 
a very fine spectral grid and then integrating the resulting spectrum. In a band model approach, 
the in-band absorption is approximated based on statistical assumptions regarding line positions 
and overlap. Temperature dependent band model parameters are computed from an atomic and 
molecular transition line atlas such as HITRAN. [See L.S. Rothman, CP. Rinsland, A. 
Goldman, S.T. Massie, D.P. Edwards, J.-M. Flaud, A. Perrin, V. Dana, J.-.Y. Mandin, J. 
Schroeder, A. McCann, R.R. Gamache, R.B. Wattson, K. Yoshino, K. Chance, K.W. Jucks, L.R. 
Brown, V. Nemtchinov, and P. Varanasi, The HITRAN Atomic and molecular Spectroscopic 
Database and HAWKS (HITRAN Atmospheric Workstation): 1996 Edition, J. Quant. Spectrosc. 
Radiat. Transfer, 60, pp. 665-710 (1998)]. These parameters define an effective single-line for 
the interval, characterized by its absorption line strength and half-width parameters, and the 
effective number of lines, n e ff. 

In MODTRAN4, the finite spectral bin single-line Voigt equivalent width W si is 
computed to determine the absorption of the effective average line. It is calculated as the 
difference between the total equivalent width, computed using the Rodgers-Williams formula 
[see: CD. Rodgers and A.P. Williams, "Integrated absorption of a spectral line with the Voigt 
profile", J. Quant. Spectrosc. Radiat. Transfer, 14, pp. 319-323, 1974], and the absorption due to 
the two line tails falling outside of the spectral band. The line tail calculations are performed for 
a line centered 0.2 cm" 1 from the edge of the 1.0 cm" 1 spectral bin; offsetting the location of the 
effective line from the center of the bin gives a more representative result for the average 
absorption of the two line tails. MODTRAN4 computes the line tail absorption by modeling the 
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tail line-shape as being inversely proportional to the square of the line center displacement, i.e., 
oc (A v)" 2 . With lines centered 0.2 cm" 1 from the edge of the spectral band, Doppler contributions 
to the line tails are small and the Lorentz line-shape denominator is dominated by the line center 



MODTRAN4 1.0 cm" 1 band model line tail absorption is defined as the absorption from 
molecular transitions centered outside of the 1.0 cm" 1 band but no more than 25 cm" 1 from band 
center. Contributions from beyond 25 cm" 1 are only considered for H 2 0 and C0 2 , and modeled 
as continua (based on the approach in FASCODE). The line tail absorption is calculated from a 
database of temperature dependent 0.25 cm" 1 integrated Lorentzian absorption coefficients. The 
line tail spectral dependence is assumed to be relatively flat so that the absorption coefficients 
can be modeled as constant over the 0.25 cm" 1 spectral grid. To justify this assumption and 
enable line tails to be modeled as Lorentzian, atomic and molecular transitions centered too close 
to a 1.0 cm" 1 spectral band edge are translated inward. A small line-shift correction is applied in- 
band to preserve the total integrated line strength. The overall error introduced into the 1.0 cm" 1 
band model by shifting line centers is small. 

MODTRAN4 computes the total 1 .0 cm" 1 spectral band transmittance T v for the n eff 
identical lines by assuming line overlap characteristic of randomly distributed lines within a 
spectral interval. Plass [see: G.N. Plass, "Models for Spectral Band Absorption", J. Opt. Soc. 
Am,, 48, pp. 690-703, 1958] showed that the transmittance due to randomly distributed identical 
lines is given by the expression 



displacement term for Lorentz half-widths less than about 0.1 cm" 1 . 




(1) 
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The Plass transmittance reduces to exact expressions in the limit of a single line, 
T v (n ef f=l) = 1-JF,//Av, and in the many line limit, T v (n efr ->cc) = exp ( - n eff W s i/Av). 

As spectral resolution narrows, direct application of the MODTRAN4 band model 
becomes inaccurate. With a A v= 0.1 cm" 1 spectral bandwidth, for example, and an effective 
5 average line positioned 0.25 A v from the bin edge, line tail absorption at 1 atm pressure contains 
significant Voigt contributions, and the MODTRAN algorithm is not applicable in this regime. 

SUMMARY OF THE INVENTION 

SERTRAN, a Spectrally Enhanced Resolution extension of MODTRAN, has been 
10 ?|| developed and its radiative transfer algorithm is described and claimed herein. SERTRAN 
W retains all MODTRAN4 capabilities while providing a factor of up to ten improvement in 
Hf spectral resolution. Due to its finer spectral resolution, SERTRAN improves treatment of line 
! s " correlation and increases overall accuracy when compared to MODTRAN4. Because of the 

J increased computational complexity, SERTRAN runs somewhat slower than the MODTRAN4 1 

n! 

1 5 CS cm* 1 binned band model, but is still considerably faster than LBL radiative transport models. 
N The MODTRAN4 and SERTRAN band models do share major radiative transport 

elements. The basic quantities computed by the band models are individual species spectral 
transmittances through homogeneous path segments. These segments are defined within single 
atmospheric layers if the Correlated-A: (CK) [see: L.S. Bernstein, A. Berk, D.C. Robertson, P.K. 

20 Acharya, G.P. Anderson, and J.H. Chetwynd, "Addition of a Correlated-*; Capability to 
MODTRAN," Proceedings of the 1996 Meeting of the IRIS Specialty Group on Targets, 
Backgrounds, and Discrimination, Vol. Ill, pp. 249-258, 1996] algorithm is selected and for 
Curtis-Godson averaged paths [see: A.R. Curtis, Contribution to a Discussion of "A Statistical 
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Model for Water Vapor Absorption," by R. M. Goody, Quart. J. Roy. Meteorol. Soc. 78, pp. 
638-640, 1952, and see: W.L. Godson, "The Evaluation of Infrared-Radiative Fluxes Due to 
Atmospheric Water Vapor," Quart. J. Roy. Meteorol. Soc. 19, pp. 367-379, 1953] otherwise. 
The individual species transmittances are themselves computed as a product of three terms: the 
5 molecular continuum (Av > 25 cm" 1 ), line tails and line centers. Only the line tail and line center 
calculations differ between the MODTRAN4 and SERTRAN band models. 

This invention features a band model method for computing individual atomic and 
molecular species spectral transmittances through a gaseous medium, comprising: dividing the 
spectral region being considered into a number of spectral bins, each having a width of less than 
10 »II 1.0 cm -1 ; calculating the equivalent width of atomic and molecular transitions centered within 
W each spectral bin; and calculating line tail absorption within each spectral bin from atomic and 
% molecular transitions not centered within the bin. 

r The spectral bins may have a width of about 0. 1 cm" 1 . The calculating step may include 

|J an exact expansion for calculating the bin Voigt equivalent width of atomic and molecular 
15 3 transitions whose centers lie within each spectral bin. The exact expansion is preferably an exact 
M modified Bessel functions expansion. The calculating step may include subtracting the line-tail 
absorption as calculated from the column strength, the Lorentz half-width, the Doppler half- 
width, and the line tail spectral displacement. The calculating step may include determining the 
Voigt line-shape function at specific spectral frequencies. 
20 The line tail calculation step may include calculating line tail absorption within each bin 

from atomic and molecular transitions centered outside of the bin using Pade approximant 
spectral fits to Voigt absorption coefficient curves. The line tail absorption calculation step may 
include determining a database of temperature and pressure dependent Pade approximant 
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spectral fits to Voigt absorption coefficient curves. There may be five Pade parameters. The 
Pade parameters may be determined from summed line tail spectral absorption coefficients. One 
Pade parameter may be determined at the center of the bin, and one at each edge of the bin, and 
one may be the derivative of the absorption coefficient with respect to the normalized spectral 
variable at the line center, and one may be the integral of the spectral absorption coefficient over 
the spectral band. The Pade parameter database is preferably generated for a plurality of 
temperatures and pressures. 

Another embodiment of the invention features a band model method for determining the 
contribution of line centers to the computation of individual atomic and molecular species 
spectral transmittances through a gaseous medium, comprising: dividing the spectrum being 
measured into a number of spectral bins; and calculating the bin Voigt equivalent width of 
atomic and molecular transitions centered within each spectral bin from an exact expansion. 

Yet another embodiment features a method for determining the contribution of line tails 
to the computation of individual atomic and molecular species spectral transmittances through a 
gaseous medium, comprising: dividing the spectral region being considered into a number of 
spectral bins; and calculating line tail absorption within each bin from atomic and molecular 
transitions centered outside of the bin using Pade approximant spectral fits to Voigt absorption 
coefficient curves. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Other objects, features and advantages will occur to those skilled in the art from the 
following description of the preferred embodiments and the accompanying drawings in which: 
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Figure 1 displays molecular transition spectral absorption near the ground: (a) Spectral features 
within a 1.0 cm" 1 bin; (b) Spectral features within a 0.1 cm" 1 bin. 

Figure 2 is a comparison of Band Model and Line-By-Line Approaches. LBL models calculate 
high spectral resolution absorption (solid curve). Band models statistically determine 
5 the average integrated strength (dash line) based on line strength, line density, and 

line width parameters. 

Figure 3 shows a High Spectral Resolution Line Tail Absorption Curve (thick dash) and Pade 
Approximant Fit (thin line) for Lorentzian H 2 0 Lines at 305K and latm Pressure. 
W The spectral curves from largest contributors are also illustrated (lines with diamond 

10 ul for lower frequency transitions and lines with squares for higher frequency 

transitions). 



§.i Figure 4 is a High Spectral Resolution Line Tail Absorption Curve (dash) and Pade 

; P Approximant Fit (solid) for Voigt lines with Strong Doppler Influence. The lower 

ity. 'i 
||.| 

curve is the residual between the upper curves. 



1 5 Figure 5 shows the in-band transmittance as a function of bandwidth for a pure Lorentz line 
with half-width 0.1 cm" 1 . Results are illustrated for line center optical depths of 0.04 
(squares), 0.2 (circles), 1 (up_arrow), 5 (down_arrow), 25 (diamond) and 125 
(left__arrow). The thick curves with symbols are the results from the modified Bessel 
function expansion truncated at I 6 , and the thin curves without symbols are the 

20 results computed assuming line tails are inversely proportional to frequency squared. 
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Figure 6 shows the in-band transmittance as a function of line center absorption for a pure 
Lorentz line with half-width 0.1 cm" 1 . Results are illustrated for bandwidths of 6.4 
(squares), 1.6 (circles), 0.4 (up_arrow), 0.1 (down_arrow) and 0.025 (diamond) cm" 1 , 
and for modified Bessel function expansions up to h (thin_dash), h (solidjine) and h 
5 (line_with_symbols). 

Figure 7 shows Doppler total equivalent width predictions as a function of line center 

absorption. The modified Bessel function expansions truncated at I 5 and I 7 , and a 
simple weak-line to strong-line interpolation formula is compared the exact result. 
Q Note the curves have been plotted with a break in the ordinate scale at 0.08 in order to 

1 0 W accentuate the residuals. 

81 

W Figure 8 shows Validation Results for MODTRAN (left) and SERTRAN (right). The band 

H model spectral transmittances are compared to FASE LBL predictions for a 0.5-km 

% path at the ground. The results from all three models were convolved with the same 2 

W 

r| cm spectral response function. 

15 Figure 9 is a comparison of FASE and SERTRAN 0.2 cm" 1 Spectral Transmittances near 2.88 
um. 

Figure 10 shows FASE and SERTRAN 0.2 cm" 1 Spectral Transmittances for C0 2 at 40-km 
Altitude and a 100-km Constant Pressure Path. The results are illustrated for the 
entire 4.3-um band (left) and an enlargement for the center of the band (right). 
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Figure 1 1 shows LWIR Validation of SERTRAN (left) and FASE (right) Spectral Radiances 
Against HIS Spectrometer Measurements for Nadir Viewing from 20-km Altitude 
over the Pacific Ocean (14 April 1986). SERTRAN radiances were convolved 
assuming a 0.364 cm" 1 (FWHM) Sine function. 

Figure 12 shows SWIR Validation of SERTRAN (left) and FASE (right) Spectral Radiances 
Against HIS Spectrometer Measurements for Nadir Viewing from 20-km Altitude 
over the Pacific Ocean (14 April 1986). 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Two RT algorithms have been implemented in SERTRAN. The two options are referred 
to as the "single-bin" and "double-bin" approaches. The "single-bin" approach is preferred for 
highest spectral resolution calculations, and it uses the traditional partitioning of the atomic and 
molecular line absorption contributions to a spectral bin. Line center absorption is defined as the 
absorption from all lines centered within the bin, and the line tail absorption is defined as the 
contribution from all lines centered outside of the bin, but within a pre-selected cutoff distance 
(25 cm" 1 in MODTRAN and SERTRAN). As is common in band model theory, an effective 
average line is defined for calculation of the line center absorption. The finite bin equivalent 
width of that line is computed with the line centered one-quarter of a bin-width from one bin 
edge and three-quarters of a bin-width from the other edge. New algorithms described below are 
required to accurately and efficiently compute the finite-bin equivalent width for bin-widths of 
size comparable to the line width. The temperature and pressure dependent spectral absorption 
from line tails is pre-computed and fit to Pade approximants. In the line tail computations, lines 
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centered too close to a spectral bin edge, i.e., within a quarter of a bin-width, are translated to the 
quarter bin-width position in order to avoid modeling their line centers as a line tails. 

The alternative "double-bin" band model approach was introduced into SERTRAN to 
simplify the line center calculations and improve upon the line tail algorithm. The SERTRAN 
"double-bin" approach defines the band model line center absorption as the in-bin absorption 
from all lines centered in the selected bin and all lines centered within half a bin-width of the 
selected bin. With this definition, the line center absorption of each atomic and molecular 
transition is partitioned between two spectral bins, the one in which the line is located and its 
nearest neighbor - hence, the name " double-bin" . On average, half the lines contributing to the 
interval are centered in the interval and half are centered outside the interval. Thus, the finite-bin 
equivalent width is computed for a line centered precisely on the edge of the spectral bin. For 
0.1 cm" 1 spectral bins, the total in-band absorption from the edge-centered line equals one-half 
the difference between the total equivalent width of the line and the 0.1 cm' 1 tail. This simplifies 
the line center calculation because only a single line tail must be computed, and, more 
importantly, the line tail begins a full bin-width from line center instead of at the quarter of a bin- 
width point. The other advantage of the " double-bin" over the "single-bin" approach is that the 
line tail absorption by definition only includes the contribution from lines centered more than 
half a bin-width from a given bin. Thus, the position of line centers does not have to be 
translated to avoid modeling line centers as line tails. 

The disadvantage of the "double-bin" approach is that the line center absorption from 
lines which are narrow compared to the bin-width is evenly partitioned between two bins when 
in reality the absorption essentially occurs in a single bin unless the line happens to be positioned 
at the bin edge. The consequence of this incorrect partitioning is that the resolution of the band 
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model is reduced. Thus, a 0.1 cm" 1 "double-bin" band model actually provides spectral 
resolution comparable to a 0.2 cm" 1 "single-bin" band model. For this reason, highest spectral 
resolution calculations should be performed with the "single-bin" method. Of course, the bin- 
width in the double-bin method could be reduced to 0.05 cm" 1 to be comparable to the single-bin 

5 0.1 cm" 1 implementation, but this would double the size of the required databases and eliminate 
much of the computation simplicity of the line center computation (because the line tail offset 
has been reduced by a half). 

Since atmospheric molecular line half-widths can be as large as 0.13 cm" 1 , the 0.1 cm" 1 
line tails used in the double-bin approach as implemented in SERTRAN (or worse yet the 0.025 

10 J cm" 1 line tails required for the single-bin approach) cannot be modeled as simply being inversely 

?. ; 'i 

W proportional to the square of the line center displacement. A description of the upgraded 
% algorithm for line tail absorption follows. 

;|; SERTRAN band model line tail absorption 
1 5 0 For SERTRAN, with its Av = 0. 1 cm" 1 band model, the MODTRAN4 line tail algorithm 

^ is inappropriate. Line tail spectral absorption cannot be modeled as constant or even Lorentzian 
so close to line center. SERTRAN greatly improves upon the MODTRAN approach by defining 
a database of temperature- and pressure-dependent Pade approximant spectral fits to Voigt 
absorption coefficient curves, k{S v ). Five parameters, k 0 , h, k 2 , jq and x 2 , are used in the Pade 
20 fits: 

k{d v ) k 0 +k l S v + k 2 5 V 2 . - = v - v cen (2) 
P l + Xi S v + x 2 S v 2 ' v AW2 
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In this parameterization, v cen is the central frequency of the spectral band and the range of the 
normalized spectral variable 8 V is from -1 to +1. Since the Lorentz line shape is proportional to 
pressure to leading order, the absorption coefficient over pressure P is fit to the Pade 
approximant. This form for the tail contributions yields an exact spectral fit in the limit of a 
5 single Lorentz tail. The five Pade parameters are determined from the summed line tail spectral 
absorption coefficients at the center and edges of the spectral band, &(0), k{-\) and k(+l), the 
derivative of the absorption coefficient with respect to <^at v cem k'(0), and the integral of the 
spectral absorption coefficient over the spectral band, k(8^ d5 v . The values of x\ and X2 
are restricted to insure no singularities for S v in No exact solution exists for the Pade 

10 ]zi parameters in a very limited number of cases, approximately 1 in 100,000. Essentially, the 
m spectral absorption coefficient data are too restrictive to allow the integral condition to be 
Q satisfied. In these cases, the derivative constraint is relaxed to enable the integrated absorption 
« coefficient to be preserved. 

4; Generally, the Pade spectral fits are extremely accurate. An example is shown in Figure 

15 J:;j 3. In this case, hundreds of molecular transitions contribute to the 305K and 1 atm pressure H 2 0 
^ line tail absorption coefficient spectrum between 200.0 and 200.1 cm" 1 , the thick dashed curve. 
The six largest contributors to the curve from lines centered above 200.15 cm" 1 are shown as 
lines with squares, and the two largest contributors from lines centered below 199.95 cm" 1 are 
shown as lines with diamonds (note the break and change of scale in the ordinate axis). As 
20 described above, the Pade fit fixes the spectral curve values at 200.00, 200.05, and 200.10 cm' 1 , 
the spectral derivative at 200.05 cm" 1 , and the integrated value. In this case, as with most, the 
spectral fit lies directly on top of the high-resolution spectrum. 
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The spectral fit to the line tails is least accurate but still quite acceptable when the 
Doppler contribution is strong. The Pade form is unable to completely mimic the exponential 
decay of the Doppler line shape. For SERTRAN molecular bands, the worse fits result from 
near-IR 0 2 lines and near-IR and visible H 2 0 lines. In the example of Figure 4, the Doppler 
5 half-widths for the contributing H 2 0 lines are near 0.019 cm" 1 and the atmospheric pressure is 
0.1 atm. Although the Pade fit satisfies the 5 prescribed conditions, small but apparent residuals 
result. 

In the preferred embodiment, the SERTRAN Pade parameter database is generated for 
six temperatures from 180 to 305K, and for pressures of 1.0 and 0.1 atm. In SERTRAN, the user 
10 || can set the parameter which defines the number of monochromatic line tail absorption coefficient 

W calculations performed per 0.1 cm" 1 spectral bin; it is nominally set to five. The line tail optical 

81 depth is modeled as varying linearly between spectral points. The monochromatic line tail 

O 

p ' optical depth and thermal emission calculations can be summed through multiple atmospheric 
% layers, with the in-band (0.1 cm" 1 ) transmittance integrated as the final step. 
1 5 Pj In principle, one could perform multiple scattering calculations at the spectral resolution 

U defined by the sub-division of the spectral bin; however this procedure would be computationally 
expensive and inconsistently mesh with the band model line center optical depths, which are 
already spectrally averaged. Thus, unless the correlated-* algorithm is selected, the multiple 
scattering calculations are performed with in-band optical depths defined for each layer. 

20 If a correlated-* option is invoked for computing molecular absorption, then the spectral 

line tail optical depth must be distributed among the ^-intervals. The procedure for partitioning 
the line tails is not well defined. The mapping from absorption coefficient value to frequency is 
not necessarily fully correlated for multiple path segments, and the mapping is generally lost in 
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the construction of ^-distributions. In fact, in the MODTRAN4 correlated-A; algorithm, the k- 
distributions are statistical representations based on the band model parameters; no exact 
mapping back to frequency space even exists. However, since the band model is based on a 
random distribution of line centers, a random partitioning of the line tail absorption among the k- 
intervals is a consistent approach. In SERTRAN, the partitioning is fixed, not random, but it is 
defined to insure that the strongest line center absorption features are well distributed over the 
spectral-bin. 

, Generalizing the Ladenbiirg and Reiche Equivalent Width Formula for Treatment of Voigt 
I Lines and Finite Spectral Bins 

I Introduction 

I Band model radiative transport models such as MODTRAN and SERTRAN [see A. 

Berk, L.S. Bernstein, G.P. Anderson, P.K. Acharya, D.C. Robertson, J.H. Chetwynd, and S.M. 

\ Adler-Golden, MODTRAN Cloud and Multiple Scattering Upgrades with Application to 

j AVIRIS, Remote Sens. Environ., 65, 367-375 (1998), and see A. Berk, P.K. Acharya, L.S. 

' Bernstein, G.P. Anderson, J.H. Chetwynd, and MX. Hoke, "Reformulation of the MODTRAN 
band model for higher spectral resolution," in Algorithms for Multispectral, Hyperspectral, and 
Ultraspectral Imagery VI, Sylvia S. Shen, Michael R Descour, Editors, Proceedings of SPIE Vol. 
4049, 190-198 (2000)] provide accurate and efficient methods for computing in-band 
atmospheric transmittances, radiances and fluxes. These models pre-process molecular transition 
data like that in the HITRAN database [see L.S. Rothman, CP. Rinsland, A. Goldman, S.T. 
Massie, D.P. Edwards, J.-M. Flaud, A. Perrin, V. Dana, J.-.Y. Mandin, J. Schroeder, A. McCann, 
R.R. Gamache, R.B. Wattson, K. Yoshino, K. Chance, K.W. Jucks, L.R. Brown, V. Nemtchinov, 
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and P. Varanasi, The HITRAN Molecular Spectroscopic Database and HAWKS (HITRAN 
Atmospheric Workstation): 1996 Edition, J. Quant. Spectrosc. Radiat. Transfer, 60, pp. 665-710 
(1998)] to define temperature-dependent total line strengths, temperature-dependent effective 
numbers of lines, an average collision (Lorentz) half-width at standard pressure, and line-tail 

5 absorption coefficients for each molecular absorption source in finite spectral bands. These 
parameters along with the Doppler half-width are used to define statistically averaged molecular 
absorption lines within homogeneous path segments. The path segments are specified by their 
temperature, pressure, molecular abundances and path length. The in-band absorption or 
equivalently definite bin Voigt equivalent width of the statistically averaged molecular lines is 

lOfi computed and combined with line position and strength distribution assumptions [see: G.N. 

W piass, Models for Spectral Band Absorption, J. Opt. Soc. Am., 48, 690-703, (1958), see: D.C. 

Ill 

!1 Robertson, L.S. Bernstein, R. Haimes, J. Wunderlich, and L. Vega, 5 cm 1 Band Model Option to 
f LOWTRAN 5, Appl. Opt., 20, 3218-3226 (1981), and see: R.M. Goody and Y.L. Yung, 
% Atmospheric Radiation: Theoretical Basis, 2 nd ed., (New York, Oxford University Press, 1989)] 
153 to determine the in-band transmittance. 
H The band model approach provides a major advantage over the more rigorous line-by-line 

(LBL) technique in that in-band values are computed directly. In LBL calculations, spectral 
optical depths are determined at a very high spectral resolution (typically, 0.01 to 0.0001 cm" 1 ) 
and then transmittances are spectrally integrated to obtain the in-band values. At each spectral 
20 point, the absorption from all contributing atomic and molecular lines must be computed and 

summed. In the long wave infrared, hundreds-to-thousands of lines often contribute significantly 
to a single spectral point. It is for this reason that a statistical approach is beneficial and often 
crucial for solving real world problems. 
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With the increasing signal-to-noise, spectral resolution and data flow of current and 
developmental infrared optical sensor, there is a desire and requirement to improve band models 
so that finer spectral resolution and higher accuracy predictions can be quickly generated. A 
critical aspect of this development is the calculation of accurate equivalent widths. The standard 
approach for calculating finite spectral bin equivalent widths involves computing the total line 
absorption i.e., total equivalent width, and then subtracting the two line-tail contributions. At 
moderate spectral resolution (> 1.0 cm" 1 ), the total equivalent width can be computed by 
interpolating between the pure Doppler and Lorentz total equivalent widths [see: CD. Rodgers 
and A.P. Williams, Integrated absorption of a spectral line with the Voigt profile, J. Quant 
Spectrosc. Radial Transfer, 14, 319-323, (1974)], and the line-tail contributions determined 
from the reciprocal-frequency squared drop off of the line shape. 

The development of a higher spectral resolution ( = 0.1 cm" 1 ) band model requires highly 
accurate total and line-tail equivalent width calculations. When the frequency displacement from 
line center to line-tail is small compared to the Lorentz or Doppler half widths, the magnitude of 
the two terms being differenced is comparable. Furthermore, the line shape of the tails exhibits 
strong Voigt behavior. In this section, an exact modified Bessel functions expansion for the 
Voigt total and line-tail equivalent width is derived. Equations for the expansion coefficients are 
given for the pure Lorentz limit, for the pure Doppler limit and for Voigt line-tails with zero-to- 
moderate Doppler contributions. 

The Problem 
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The objective of this section is to determine an accurate and efficient algorithm for 
evaluating the Voigt line shape absorption W L = W^(Su 9 y c9 y d ) integrated between spectral 
frequencies v' and co for a molecular transition centered at frequency v 0 : 

W A - J[l-exp(- &/„)]<* v 



(3) 



5 where the Voigt line-shape f v is the convolution of Lorentz and Doppler line-shape functions: 



f = Js ~ f 



expl 



-dt 



10 1 



In these equations 

y c > 0 is the collision (Lorentz) half-width at "half maximum [cm* 1 ], 
y d > 0 is the Doppler half-width at "1/e " of the maximum [cm' 1 ], 
A = v f - v 0 is the spectral frequency displacement from line center [cm" 1 ], 
jS{ v is the spectral frequency integration variable [cm" 1 ], 

S is the line strength [atrn 1 cm' 2 ], and 
Q u is the column density [atm cm]. 

■ttv 

m Since the Voigt line shape function f v is symmetric in v , the absorption term W_ K equals 

1 5p 2FF 0 - , and it suffices to restrict the analysis to A > 0 . For A > 0 , W A is the (single-sided) 
line-tail absorption. 



(4) 



The Ladenburg and Reiche function 

Ladenburg and Reiche [see: R. Ladenburg and F. Reiche, Uber selektive Absorption, 
20 Ann. Phys. 42, 181-209 (1913).] derived the exact expression for the total absorption (N = -qo) 
from a Lorentzian line in terms of Bessel functions of the first kind with imaginary arguments: 



W 



Lorentz 



= 2W 0 Lore '" z =Suexp 



Su 
2ny c 



( Su ^ 



f Su ^ 



\2ny c ) \l7ty c ) 



(5) 
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(The modified Bessel functions, I Q and I, , can be evaluated using the polynomial 
approximations listed in Section 9.8 of Abramowitz and Stegun) [see: M. Abramowitz and LA. 
Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical 
Tables, (New York, Dover Publications, 1965).] Ladenburg and Reiche also showed that the 
5 line-tail absorption W A orentz in the limit of A » y c , i.e., where f v =y € jKV 2 , can be expressed in 

terms of an error function. In the next few sections, an expansion is derived for W A orentz that is 
exact for any spectral displacement A . 



p Strategy for Determining Wn 
10 If In this section, the Ladenburg and Reiche formula is generalized to provide an exact 



3 I I 



m modified Bessel function expansion for Voigt line-tail absorption. The line-tail absorption, Eq. 
Pf (3), can be integrated by parts by setting P = 1 - exp(-Suf v ) and dQ = dv and the integration 
£1 variable switched from the frequency variable v to the Voigt line-shape function itself: 

■at? 

!l W A = Su f \v f exp(- Suf)df -A[l- exp(- Sfc/J] (6) 

Ci 0 
15 Defining angle 0 by the relationship f(d) = / A (l - cos#)/2 and substituting into Eq. (6) gives 

JF A = Suf A exp(^- y- A J Jv^ sin [-f ] cos (f") eXp (^ ^ A cos ^)^ (7) 

The coefficient of the exponential in the integrand is bounded for all 0 between 0 and 
7t. At 0 = nj equals f A , which implies that v 0=jr equals A and the coefficient equals zero. At 
0 = 0 , the situation is more complicated. Since the Voigt function / is zero at 9 = 0, the 
20 spectral frequency v approaches infinity. However, this is exactly compensated by the sine 

term. With the boundedness of the coefficient established, a generalization of the Ladenburg and 
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Reiche expression is obtained by expanding the coefficient of the exponential in the integrand of 
Eq. (7) (multiplied by n f h ) in a Fourier series: 

xf L v e sin^cosf = ^ + EF.(A)cosM (8) 

II I n =] 

where 

n 

5 V n = V n {A) = h \v e sm6coine)de (9) 

0 

Substituting Eq. (8) into Eq. (7) and noting the integral representation of the modified Bessel 
functions of integer order [Abramowitz and Stegun, Eq. 9.6.19], the desired expansion is 
6? obtained: 

I FK 4 =S«exp(-y/ A ) 

S % K 

10 The higher order modified Bessel functions can be computed from I 0 (z) and I^z) using the 
Q recurrence relationship [Abramowitz and Stegun, Eq. 9.6.26], although for small z , the 
Rl ascending series [Abramowitz and Stegun, Eq. 9.6.10] is numerically more stable. With the 
j3 Bessel functions defined, the calculation of Eq. (10) reduces to determining the Voigt line shape 
function at displacement frequency A , i.e., f A and to determining the Fourier coefficients V n . 

15 

Computing the Fourier Coefficients 

Determining the Fourier coefficients directly from Eq. (9) is complicated by the spectral 

frequency term. For the pure Doppler and Lorentz limits, the inversion is straightforward. 
However, in general, the calculation of the Fourier coefficients is greatly simplified by returning 
20 to spectral frequency as the dependent variable. This is accomplished by first substituting for 
angle Q\ 



f/.(fA)4n/.(f/ 1 ]l-4l-exp(-&/ 4 )] (10) 
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V n = 2 J v, cosf«^ )# ; cos^ = 1 - y- C 1 1 ) 

Next, integrate by parts with Q = v f = v and rfP n = cos(n^ )<// = cos(« 0 A )<//„ 

F B =2vP„r%2]p/v . (12) 



Av 



After expanding the P„ integrals, the Fourier coefficients for n < 7 have the following form: 

V 0 =2F 0 , (13a) 
V^V,-2F X , (13b) 
^=4^-3^+16^/3 , (13c) 
V 3 = 6^-15^ +10F 0 -16F 3 , (13d) 
F 4 =8F,-28r 2 +56r,-35F 0 +256F 4 /5 , (13e) 
V 5 = 10F 4 -45F 3 + 120F 2 -210F, + 126F 0 -512F 5 /3 , (13f) 
r 6 =12F s -66F 4 +220r 3 -495F 2 +792r i -462F 0 +4096F 6 /7 , and (13g) 
F 7 = 14F 6 -91K 5 +364F 4 -lOOlF, +2002F 2 -3003F, + 1716F 0 -2048F 7 . (13h) 
where 

F„ S F„(A) S A/ A +-^]/rVv (14) 
/a a 

In this form, calculation of the Fourier coefficients reduces to determining Voigt function 
moment integrals, F„ . 



Weak Line Limit 

In the weak line limit, i.e. when Su is small, line absorption is known to be proportional 
to Su. Substituting lXSuAn) = [(SufJ4) n /n\\{l + c{(Su) 2 }} into Eq. (10), the line tail 

{H:\pa\corpV16849\00004\A0455157DOC} 23 



absorption expression reduces to the correct weak line limit. The term o(z* ) is used here and 
throughout to indicate terms of order z n and higher. In the total line absorption limit, F 0 (o) is 
one and = 2W 0 equals Su to leading order. 



5 Strong Line Limit 

When the argument of the modified Bessel function is large, the value of the function is 

independent of n to leading order in (l/z) [Abramowitz and Stegun, Eq. 9.7.1] with 
e~ z I n (z) ~ (2nz)~ l/2 . Thus, in the asymptotic limit the modified Bessel functions can be factored 
0 out of Eq. (10), leaving the sum of Fourier coefficients. Taking the0 approaches zero limit of 
10 Eq. (8), the Fourier coefficient sum becomes 

F ° (A) +i^(A) = V^A 05) 



5 ~.r 



M 2 
IU Substituting this relationship into Eq. (10) results in the square root curve of growth: 



2 ^ 



; (r c >o) (16) 



.SufJ 

For numerical implementations, it is recommended that the strong line limit be forced by setting 
15 a last Fourier coefficient by the constraint of Eq. (15). 

The Lorentz Limit 

In the Lorentz limit (y d = 0;y c > 0), the line-shape function has a simple form: 

r Lorentz Y cfa /I 

20 It is convenient to define a as the ratio of the spectral displacement from line center A to the 
Lorentz half-width y c : 
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a = A/y c 

The first six Fourier coefficients have the following form: 



(18) 



TtV, 



Lorentz 



a 



cot a + — y 



2 "' " a' +1 
n V x Lorentl = a - (a 2 - 1) cot' 1 a 



71 F, 



Lorentz 



2a 



acot~ a-l + 



1/3 



a 2 +l 



Lorentz 



a 



5a 2 +--a(5a 2 +3)cot _1 a 



Lorentz 



2a 



' i \ 4 17 2 11 1/15 

a(la 4 +Sa 2 + 2jcot _1 a - la* - — a 2 - — + 



15 a 1 +1 



Lorentz 



= a 



42a 6 + 56a 4 +^« 2 + ^-«(42a 6 +70a 4 +35« 2 +5) cot" 1 a 



(19a) 
(19b) 

(19c) 
(19d) 
(19e) 
(19f) 



For « > 2 , each of the y^""" 2 are proportional to a 3 and are small for A«y c . However, for 
larger these expressions are numerically unstable. This limit is more easily evaluated using 
asymptotic series. 

As suggested above, convergence to the strong line limit is forced by setting a final 
coefficient based on the normalization condition, Eq. (15). This is quite an important condition. 
The general approach for calculating the absorption integrals involves expanding the coefficient 
of the exponential in Eq. (7) in a Fourier series. However, in the strong line limit, the integrand 
is completely dominated by values near 9 = 0. The Fourier series fits the entire domain of 0 
values, giving no special weighting to the only part of the domain that significantly contributes 
for strong lines. The normalization constraint helps tailor the fit in this region. 
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In Figure 5, the in-band transmittance is plotted as a function of band half-width. The In- 
band transmittance increases as bandwidth increases. If the bandwidth is a factor of 3 to 4 times 
larger than the Lorentz half- width, the line tails can be reasonably well represented by the 
reciprocal-frequency squared relationship (dashed curves). For narrower band models, this 
relationship breaks down and the exact expansion is required. 

The convergence of the Lorentz equivalent width expansion is illustrated by Figure 6. In 
each of the curves, the final Lorentz Fourier coefficient is computed using the normalization sum 
Eq. (15) rather than the exact value from Eq. (19) to insure proper convergence in the strong line 
limit. In-band transmittances above 5% are accurately computed using just 3 terms in the 
! modified Bessel function expansion, up to I 2 . Adding two more terms, the expansion is 
5 accurate for transmittances above about !/ 2 %. Using all 7 terms defined above, the predictions 
! are accurate for all transmittances above 0.001. 

In generating the 7-term curves, the predicted transmittances were compared to the 

j maximum exp(- Su f£ orentz ) and minimum exp(- Su f Q Lorentz ) possible transmittance values, and 
! replaced by the average of the extreme values when the prediction lied outside this range. This 

replacement was only necessary for transmittances below 0.001 (not for values shown in Fig. 6). 

Even though convergence to the strong line limit is guaranteed by the normalization condition, 

there is an intermediate region for very small transmittances where predictions diverge slightly. 

Transmittance errors as large as 0.004 can result for transmittances near or below 0.00001 . 

Pure Doppler Limit 

In this section, the expression for the Voigt line-tail absorption, Eq. (10), is applied to a 
pure Doppler line {j c = 0). In this limit, the line shape function has the following form: 
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J V 



Doppler = exp(-fc„) . 



(20) 



Inserting this definition into Eq. (14), the Voigt function moment integrals have a relatively 
simple form: 



*p Doppler _ 



expl 



U) 



1 + 



1 n 



2\nb, 



expl 



(nb A )erfc^jnbl 



(21) 



5 The Doppler Fourier coefficients are computed by directly substituting these moment integrals 
into Eq. (13). 

In the limit of total Doppler line absorption, the complementary error function terms are 
all one (since b 0 = 0 ), and the following expression results for the total Doppler equivalent 
W width: 



jyr Doppler _ j^Doppler 



10 u 



- Su exp 



' -Su ^ 



' Su ^ 

+ 0.58578644 7, 
-0.25271345/ 
-0.09628129/, 
-0.05165314/ 

All the coefficients in this expansion can be represented analytically. 

Applying the strong line limit, Eq. (15), to pure Doppler lines indicates that the 
y Doppier coe fg c i en t s mus t sum to zero to avoid a square root curve of growth. The coefficients 




-0.57765281/. 



-0.14605219/ 



-0.06867590/, 




(22) 
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actually included in Eq. (22) sum to -0.107. Since the sum must approach zero, terms beyond 
those listed in Eq. (22) cannot all negative. 

The accuracy of the truncated series Eq. (22) for calculation of the Doppler total 
equivalent width is illustrated in Figure 7. For comparison, results are also presented for an 
often-used simple interpolation formula between the weak and strong line Doppler limits: 



W D : Ppler * Y d ^21n(^ + i(S M /rj] (23) 

If line center absorption Su/y d ^n is less than 10, the result from using the first 6 terms in Eq. 
(22), i.e., up to I 5 , is essentially exact. For even stronger lines, errors quickly become 
appreciable. For a line center extinction of 100, the 6 and 8 term expansions produce 24.2% and 
12.5% errors, respectively. In this limit, the asymptotic expansion [see: C.B. Ludwig, W. 
Malkmus, J.E. Reardon, and J.A.L. Thomson, "Handbook of Infrared Radiation from 
Combustion Gases, Ed. R. Goulard and J.A.L. Thomson, NASA SP-3080 (1973), Eq. 2-B-2] or 
even Eq. (23) is preferable. 



Voigt Line Tails with Moderate Doppler Dependence 

The objective of this section is to obtain an analytic expression for the Fourier coefficients 

V„ of Eq. (10) for Voigt lines with at most moderate Doppler dependence. To be more specific, 

the analysis is restricted to line tails W A for which x A is large, where 

«-4 ; X ^a + b v = . (24) 

yI Yd 

From Eq. (13), the calculation of the Fourier coefficients reduces to calculation of the Voigt 
function moment integrals of Eq. (14), F„. To derive an expression for these integrals, a 
workable form for the Voigt line shape must be derived. It follows from Eq. (4) that the Voigt 
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line shape is proportional to the real part of the complex probability function w(z v ) [Abramowitz 
and Stegun, Eq. 7.4.13]. From the continued fraction representation for w(z v ) [Abramowitz and 
Stegun, Eqs. 7.1.4 and 7.1.15], the Voigt integral can be expanded in an asymptotic series: 



ajx v 



We 



1 + - 



d. 2 ± 

2 x. 



1 



2\ 



4 x„ x 2 



v J 



1^-105-^ + 210^-120^ + 

8 x v x v x ) 



ok?) 



(25) 



The expression for f v can now be raised to the n th power, expanded in powers of \jx v , and 
integrated to give the result: 



- J 

arc tan (a) 



COS 2 "0> - 



n(n + 1) 



2(n + 2) 



2\ 



l + a 



___2«+2 l 

COS f 



+ 



n{n + \\n + 2\n + i) 



l + a 2 



, v , cos 2 "^ 
8(n + 3X« + 4) ^ x A ) 

n{n + ljn* + 26n 3 + 251k 2 + 874n + 1008) f 1 + a 
48(n + 4X" + 5X« + 6) ( x A 



cos <p 



d<j> 



(n + l)a/x A ( n + 6 n 2 +8n+l9 ) (n + l)a/x 



(« + 2Xi+« 2 r i i+a2 2 (" +3 ) j(«+4Xi+« 2 r 

^ 2(» 2 +17n + 90) _ 5« 3 +99n 2 +706« + 1800 
s(l + a 2 } 6(» + 5Xl + « 2 ) 

7n 4 + 144» 3 + 1 127w 2 4- 3942n + 5400 
h 24(n+4X« + 5Xl + a 2 j 



(n + eil + a 2 } 



(26) 



where a is the ratio defined in Eq. (18) and the cosine integrals can be analytically evaluated: 



x/2 

Jcos 

arctan(a) 



2m 



(2m)! f a Si (r!) 2 ( 4 Y 

"4 m M 2 r 0t a \ + a 2 k{2r + ^K\ + a 2 ) 



3j(2r + l)! 



(27) 
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These expressions can now be substituted into Eqs. (13) and (14) to compute the Fourier 
coefficients. 

If the Lorentz half-width is large compared to the Doppler half-width, then x A is large even 
when A is zero and the equations of this section can be used to evaluate the total Voigt 
5 equivalent width. In this limit, the expression for the Voigt function moment integrals simplify 
considerably and the Voigt equivalent width series to order 0{x 4 ) actually truncates after a finite 
number of terms. 

Summarizing the Algorithm 
10 J A generalization of the Ladenburg and Reiche absorption formula to Voigt lines and their 

W line-tails has been derived. Given the column strength Su , the Lorentz half-width , the 

Si 
fll 

jl Doppler half-width y d and the line-tail spectral displacement A , an exact expansion for the line- 

^ tail absorption W A is given by Eq. (10). The Voigt line-shape function / A at spectral frequency 
S* A can either be computed using the asymptotic series of Eq. (25), or more accurately using the 
15 pi Humlicek rational approximation [see: J. Humlicek, "An Efficient Method for Evaluation of the 
Complex Probability Function: The Voigt Function and it Derivatives", J. Quant. Spectrosc. 
Radiat. Transfer 21, 309-3 13 (1979)]. For the special case of spectral frequency A = 0 , the 
Voigt line-shape is conveniently computed from its complementary error function representation. 

The modified Bessel function terms e~ z I n (z) for real z can be computed using polynomial 
20 approximations [Abramowitz and Stegun, Sec. 9.8], recurrence relations, and/or ascending 

series. To define the Voigt Fourier coefficients V n , the variables a, a andx v [Eqs. (18) and (24)] 
are introduced. The Voigt Fourier coefficients are defined as an asymptotic series in \jx v [Eqs. 
(13) - (14)], which themselves depend on the Voigt moment integrals, Eqs. (26) - (27). For the 
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limiting cases of pure Lorentz and pure Doppler lines, exact expressions are derived for the 
Fourier coefficients, Eqs. (19) for the Lorentz case, and Eq. (21) for the Doppler case. 

Many radiative transfer problems require calculation of the single-line fmite-spectral-bin 
equivalent width W sl . For a finite spectral interval containing a line-centered A wavenumbers 
5 from the left edge and A' wavenumbers from the right edge, the single-line equivalent width 
equals the sum of the contributions from the two sides of the line, (W Q - W A )+ (W 0 + W A ) . If the 
Lorentz half-width is significantly larger than the Doppler half-width, x 0 is large and all of the 
absorption terms can be computed using the formalism described here. On the other hand, if 
CJ x A (or x A .) is large but x 0 is not, the total Voigt equivalent width 2W Q can be estimated using 
10 ?! the Rodgers and Williams interpolation formula [C. D. Rodgers and A. P. Williams, (1974)], 
|f with the Lorentz total equivalent width computed from the Ladenburg and Reiche formula and 
52 the Doppler equivalent width obtained from either Eq. (22) above or from strong line limit 
Q expression (e.g., Su > 22.346 y d Vn) [Eq. (2-B-2) in C.B. Ludwig, et al, 1973]. Finally, if 
p! x A (or x A .) is not large, then the spectral bin half width A is comparable to or smaller than the 
15 U Doppler half-width y d . In this limit, the Voigt line-shape is relatively flat within the interval and 
a numerical integration using the Humlicek rational approximation for the line shape should 
converge rapidly with only a few evaluations. 



SERTRAN Validation 

20 SERTRAN validation efforts have concentrated on the LWIR through SWIR spectral 

regions, and all validations shown here were performed with the "double-bin" band model 
approach. SERTRAN was validated by comparison to predictions of the FASE LBL model and 
to measurements from the High-resolution Interferometer Sounder (HIS) Fourier Transform 
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Spectrometer [see H.R Revercomb, H. Bujis, H.B. Howell, R.O. Knuteson, D.D. LaPorte, W.L. 
Smith, L.A. Sromovsky, and H.W. Woolf,: "Radiometric calibration of IR interferometers: 
Experience from the High-Resolution Interferometer Sounder (HIS) aircraft instrument " RSRM 
'87, Advances in Remote Sensing Retrieval Methods. A. Deepak, H. Fleming, and J. Theon, Eds. 
5 1989], The FASE model has itself been extensively validated against measured data, and it 

provides the benchmark for determining the accuracy of the SERTRAN band model results. The 
HIS spectrometer measures high quality, sub-wavenumber calibrated spectral radiances between 
600 and 2700 cm" 1 (3.7 to 16.7 jim) for direct validation. 

10 |i SERTRAN to FASE validations 

SERTRAN and FASE spectral transmittance predictions were compared for frequencies 

2| between 500 and 4000 cm" 1 (2.5 to 20 |iim) and for horizontal (constant pressure) paths between 

® 0 and 40 km altitude. Path lengths were selected to produce band average transmittances near 

;f I 50% because transmittance residuals tend to be largest for transmittances near 0.5. 

15 m SERTRAN models line correlation and overlap more accurately than MODTRAN 

because of its higher spectral resolution. Also, the SERTRAN line tail algorithm improves upon 
the MODTRAN algorithm. Thus, the SERTRAN to FASE residuals will be smaller than 
MODTRAN to FASE residuals when all three calculations are degraded to a common spectral 
resolution. This is illustrated in Figure 8. Spectral transmittances are compared for the 0.5 km 

20 horizontal path at 1 atm pressure containing U.S. Standard profile H 2 0, C0 2 , and 0 3 

concentrations. High spectral resolution SERTRAN and FASE predictions were convolved with 
a non-overlapping 1 cm" 1 rectangular slit, and subsequently degraded using a 2 cm" 1 triangular 
slit to mimic the MODTRAN 2 cm' 1 result. For both MODTRAN and SERTRAN, no strong 
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bias is observed, i.e., the residual curves straddle the zero residual line. However, sporadic 
MODTRAN residuals approach 0.1 while all the 2 cm" 1 SERTRAN residuals are under 0.02. 

SERTRAN is compared to FASE at 0.2 cm" 1 spectral resolution in Figure 9. The 
scenario is identical to that of Figure 8, a /2-km ground path with H 2 0, C0 2 , and 0 3 . The 
5 spectral range has been reduced so that spectral details can be resolved. For much of the 
spectrum the residuals are under 0.01, but occasional spikes are large as 0.07 are evident, 
consistent with the accuracy of the equivalent width residuals calculations. 

Basic to the SERTRAN band model is the premise that molecular line absorption is 

CI substantial within two 0.1 cm" 1 spectral bins. At higher altitudes, this assumption breaks down - 

fefl 

10 W molecular line widths are small compared to the 0.1 cm' 1 bandwidth and most of the individual 
Sf line absorption occurs in a single bin. As a result, SERTRAN essentially becomes a 0.2 cm" 1 
2 band model. This is illustrated in Figure 10. SERTRAN 0.2 cm" 1 spectral transmittances are 
Q compared to FASE results for a horizontal path at 40-km altitude and with a 100-km range. The 
Ill residuals are larger than at lower altitudes. Individual narrow yet strong absorption lines are 

15 P equally partitioned into two 0.1 cm" 1 spectral bins yielding too much absorption in one bin and 
too little in its neighbor. 

SERTRAN to HIS validations 

SERTRAN predictions were compared to HIS spectrometer measurements [see H.E. 
20 Revercomb, H. Bujis, H.B. Howell, R.O. Knuteson, D.D. LaPorte, W.L. Smith, L.A. Sromovsky, 
and H.W. Woolf,: "Radiometric calibration of IR interferometers: Experience from the High- 
Resolution Interferometer Sounder (HIS) aircraft instrument," RSRM '87, Advances in Remote 
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Sensing Retrieval Methods. A. Deepak, H. Fleming, and J. Theon, Eds. 1989] from two separate 
campaigns: an airborne nadir measurement from 20-km altitude over the Pacific Ocean (14 
April 1986) and a ground-based zenith measurement from the GAPEX experiment in Denver (31 
October 1988). For both campaigns, atmospheric temperature and H2O profiles were derived 
5 from radiosonde data. The validation results from the two sets of comparisons are similar. The 
results are presented from the nadir measurement comparison, for which a FASE validation was 
previously performed. [See: J. Wang, G.P. Anderson, H.E. Revercomb, and R.O. Knuteson, 
"Validation of FASCOD3 and MODTRAN3: Comparison of Model Calculations with Ground- 
Based and Airborne Interferometer Observations Under Clear-Sky Conditions," AppL Optics, 35, 

O 

10 jj pp. 6028-6040, 1996], 

11} 

E HIS measurements for the 600-1080 cm" 1 spectral band and the SERTRAN and FASE 

m 

S J predictions are illustrated in Figure 1 1 . The upper curves overlay the measured and modeled 
a spectral radiances; the lower curves contain the measurement minus model residuals plotted with 
M matched vertical scales. The HIS minus SERTRAN residuals are at the few percent level with 
15 Gl spectral radiances and residuals near 80 and 3 mW m' 2 sf 1 / cm" 1 , respectively. The unapodized 
p; HIS Fourier transform spectra were reported to have 0.364 cm" spectral resolution in this band. 
For these comparisons, the SERTRAN results were computed using a 0.364 cm" 1 full width at 
half maximum (FWHM) sine function. The results seem to indicate a mismatch in the filter 
function - the residuals contain more jitter than the earlier SERTRAN to FASE comparisons 
20 would suggest. Nevertheless, the agreement is good. Both SERTRAN and FASE modeled the 
sea surface as a blackbody (unit emissivity) which accounts for the baseline offset between 800 
and 1000 cm" 1 . The residual spectrum between 1000 and 1080 cm" 1 indicates that the O3 column 
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amount was overestimated in the models - the main effect of the ozone is to attenuate the surface 
signal, and lowering the column amount should eliminate the 9,6 jam O3 band spectral feature. 

In the short-wave IR, the HIS spectrometer has a band extending from 2000 to 2700 cm" 1 , 
3.7 to 5.0 |um, with a reported spectral resolution of 0.637 cm" 1 . Below around 4.0 pm, the solar 
5 components become important for daytime measurements. FASE includes the solar reflection off 
the surface in its nadir atmospheric radiance calculations (ocean surface albedo set to 5% in the 
S WIR model calculations), but atmosphere scattered solar radiance contributions are neglected. 
SERTRAN, like MODTRAN, includes the atmospheric scattered contributions in its 
Pi calculations. This is evident in the SWIR HIS comparison, Figure 12. Between 2400 and 2600 

10 gl cm FASE slightly under-predicts the spectral radiance. Although the SERTRAN predictions 

W 

Si are noisier, the small bias in this region is decreased. This result had previously been 

M demonstrated using MODTRAN, albeit with the spectral data degraded to 2 cm" . The 

z m SERTRAN comparison has more jitter than the FASE comparison, especially between 2000 and 

■Is I 

jli 2200 cm" 1 , but again this is most likely due to a spectral filter mismatch. Otherwise, the 
15 £| SERTRAN and FASE results are comparable. 

The reformulation of the MODTRAN band model and development of SERTRAN together 
provide a basis for future, rapid sub-wavenumber radiative transport analysis in the terrestrial 
atmosphere. The low-altitude IR validations performed to date demonstrate that SERTRAN 
achieves MODTRAN accuracy but at higher spectral resolution, with band average 
20 transmittance residuals of order 0.01 or better. Furthermore, comparisons of SERTRAN to 

MODTRAN at a common spectral resolution (> 2 cm" 1 ) show approximately a 4-fold decrease in 
residuals. 
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